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ABSTRACT 

An infinite number of perturbed fc = Friedmann cosmologies filled with dust and radiation 
is found. As we go up the sequence, the solutions contain an increasing number of integration 
functions. With the coordinate gauge adopted to co-move with the perturbed matter, the 
solution of a pair of coupled equations for the trace h of the metric perturbations and for the 
radiation density perturbations Spr is the key to the problem. An estimate is given of the 
temperature variation of the cosmic microwave background radiation due to the Sachs- Wolfe 
effect. It is found that the temperature fluctuations can grow faster than in the absence of 
radiation. 
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Introduction 



The first-order perturbations of a Friedmann universe with a flat 3-space filled either with dust or 



radiation have been obtained by Sachs and Wolfe 1967. They computed the temperature fluctuations 
of the cosmic background radiation in the dust-filled universe, assuming that the domain in which the 
photon travels is matter-dominated. This classic prediction for the temperature fluctuations overestimates 



the experimental value ( Mather 1992 ) by several orders of magnitude. The prediction has been made by 
considering only the gravitational perturbations along the path of the photon. Obviously, the inclusion of 
random temperature fluctuations on the surface of the last scattering can only increase the effect. Neither 
did the discovery of large-scale structures (voids and walls) bring us closer to the resolution of this paradox. 



Though it has been put forward (Mather 1992) that the observed fluctuations are primordial (as opposed 
to propagation effects), this is hard to take seriously unless we are able to reduce the magnitude of the 
Sachs- Wolfe effect. 

Here we compute the perturbations of the fc = Friedmann universe in the presence of both dust 



and radiation. Perturbations of two-fluid cosmologies have been studied in the past decades (Kodama 



1987 , Koranda 1994 , Perjes and Komarik 1995 ). Instabilities are known to exist in two- fluid cosmologies 
( Mukhanov 1992 ). One could reduce the magnitude of the metric fluctuations affecting the photon orbits 
by assuming, for example, that the observed large structures are late emerging manifestations of the 
instabilities. The main result of the present paper, in Sec. 6, is an infinite series of solutions, with an 
increasing number of integration functions. Like in the pure dust model, the solutions contain a relatively 
increasing and a relatively decreasing mode of the radiation density perturbation. The growth is, however 
faster in the presence of radiation and even faster as we go up the sequence of our solutions. In Sec. 7, we 
estimate the magnitude of the Sachs- Wolfe effect on the temperature variation of the cosmic background 
radiation. We find that the relatively growing mode is incompatible with the observations already with the 
first solution. 



The energy-momentum tensor is a sum of those of the two media, 



^ b — ^ m b'-'-rb 



(1) 



The contribution of the dust has the form 



7^™fc = P™"""b- (2) 

For the radiation, 

T.-^b = IprU'^Ub - lprS%. (3) 

The four-velocities are normalized M°Ua = 1- We assume that, after the decouphng, the conservation law 
Tj°^ = (p + p) Ua-fiU^ ~ p ,a ~\-UaU^P,b = applies separately to both the matter and the radiation components: 



T° = n 



T^hn — 0. 



(4) 



This is justifiable because the decoupling occurs near the time of equal matter-and-radiation density, and 
after that, coupling occurs only via the universal gravitational interaction. We get two energy conservation 
equations by transvecting with u^. 



0. 



and two momentum conservation laws: 



i[a,6]"m = 0, 4pr {Ura.b - Urb,a) = Pr,a - Urau\.pr,b- 



(5) 
(6) 
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2. The perturbed model 

We adopt the conformal form of the metric 

gab = a^iri) {rjab + hab) , (7) 

where {rjab) — diag and hab is the metric perturbation. The indices of the perturbed 

quantities are lowered and raised by the Minkowski metric rjab — We shall not use the explicit form 
of the scaling function a(ry) as long as possible. In the perturbed universe, hab 7^ 0, the densities of the 
components can be written to first order 



(8) 



where i stands either for m (matter) or r (radiation). Here pr and pm are the unperturbed densities 



Pm{ll) = PrM-% , Pr{ll) = PrQ-^ (9) 

and 5 Pi are the first-order density perturbations. According to the unperturbed models, the radiation 
density pr dies out faster than the matter density, p„i • 

The dipole effect of the cosmic background radiation provides an experimental value for the local 
relative velocity of the order of 600 km/ sec. Thus we have good reason to assume that 5u\-^-^ ^ 5u].. Here 
we stick to a comoving gauge. For gauge invariant methods, c/. Refs. Magueijo 1993, Russ 1993, Kodama| 



1987 



Dunsby 1991, Mukhanov 1992. We choose coordinates comoving with the matter: 



where the coincident unperturbed velocities are 

a ^ xa 

The normalization conditions imply that /iqo — and Su'^ = 



(10) 
(11) 

(12) 



0. 



The permissible coordinate transformations in the comoving gauge are (3achs and Wolfe 1967) 
where the first-order function has the form 



a 



r = c"(x^) 



and a:" are the space coordinates for a = 1, 2 and 3. The metric perturbations transform 





hap 

hop 



a 



- Ca,0 

hop, 



Cp.a -+ 

hoQ 



2bllaP 



(13) 



(14) 



,0 



The velocity perturbation 



where a prime (') denotes derivative with respect to the time coordinate rj 
is gauge invariant to the required order. The density perturbations transform as follows (Brauer 1990) 



= 6pi + -p'^. 
a 



(15) 
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2.1. Energy conservation 

The energy conservation (|^) for the first-order perturbations of the dust reads 



— + 7^h\ =0 (16) 

where = ft," is the trace of the metric perturbation and a prime denotes partial derivative with respect to 
the conformal time rj. Hence we get 

5p.n = Pm [Eixl") - \h) . (17) 

Here the integration function E{x^) depends only on the space coordinates x°' [a — 1,2 or 3). A gauge 



transformation ( |13[ ) alters E as follows, 

E{x^) = E{x^)+c'^^. (18) 
The energy conservation law (|^) for the radiation has the form 

( pl V9K ) ^ = 0. (19) 



Collecting the first-order terms, we have 

'3 6pr 1 , V 



(20) 



2.2. Momentum conservation 

For the dust, the conservation law (^ simplifies, 

{ahoj = 0. (21) 

This has the solution 

ah^o = F^{x>^) . (22) 



The coordinate freedom ( |13| ) makes it possible to arrange (White 1973) f " = 0, whence 

h°% - . (23) 
The remaining gauge transformations are still of the form (p^), with 

Ab = 0. (24) 

Here the Laplacian is defined by A6 = —ri"^b,af3 ■ 

The momentum conservation law (^ reads for the radiation perturbations 

.'/^5ury ^^ap-'/^6pr^^ . (25) 
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2.3. The potential v 



In (Perjes and Komarik 1995), the potential v has been introduced under the assumption that the 
vorticity of the radiation vanishes. We now show that this potential exists for generic perturbations, without 
resorting to any assumption about the vorticity. We introduce the potential v = v (r],x^) by writing 



Under gauge transformations, 

Thus we can integrate Eq. (p5|) as follows: 



jap^^^'^Spr = v'. 



(26) 



Pr 



(27) 



such that the function of integration has no significance in (^6|) and is dropped. Hence the velocity 
perturbation 5ura = 5 {gabU^) has the form 



rjapSu'^. = Pr^^'^V^a (??, x) - Fa (x) 

Substituting and Spr in the perturbed energy conservation [Eq. (pO|)] , we obtain 

3v" - Aw + yii^aoh' 0. 



(28) 



(29) 



3. Einstein equations 

The perturbed Einstein tensor will be written 

Gl = oG% + SG% . (30) 

Here oG°^ is the unperturbed tensor and SG^ the first-order part. The Einstein equations for the first-order 
quantities are 

6G\ = -iSpr+Sp,n) (31) 

SG^o = -fapr&t" (32) 
5G"^ = ^S'^^Spr. (33) 

Substitution of the ^achs and Wolfe 1967] expressions for SG"], and separating the trace-free part of the 
metric perturbation 

Sap = haf3 - IVaph (34) 

yields 

.fJ.i' ^ 3 



5^".^, + |A/i - 2^h' - - 2a^ {Spr + Spm) , (35) 



-5*"^.^ ' - l^'"' + ^'^"^ - 4(^2^-iipj/t°"==: l^^prK , (36) 
2h" + 4^/i' - 5^"^ „^ - f A/i = -2a^5pr , (37) 
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-i/i,"^-i<5"^A/i. (38) 
The sum of (BSl) and (p^) gives the simple relation 



h" + sLh' = -a^ {2&pr + 5p,n) ■ (39) 

Taking the divergence of Eq. ( p6| ) and subtracting the 77 derivative of (p5|), we get the integrability condition 
of these equations. With the help of Eqs. (p7|), ( ^9[ ) and (|2^), we may verify, however, that this integrability 
condition is satisfied identically. 



4. Plan of solution 

In this section we review the steps of the solution procedure assuming no a priori knowledge of the 
scaling function a. First, we write the radiation velocity perturbation 5ura = 5 (gabU^) as 

6ura = a^riaf3Suf^ + ahoa- (40) 

Eliminating Su" from Eqs ( p7|) and (p8|), we obtain 

{l^ + h}" -h^ = 0- (41) 

In the first place, the coupled system consisting of this and Eq. (^) is solved for the radiation density 
perturbation Spr and for the trace h. We next express (5uf from Eqs. ( ]25| ) and (pO[): 

77Q/3(5uf = ifl^^Pr J ap^^^^'^Spr^adri -a^^hoa- (42) 



We determine the divergences of the trace-free part from Eqs. (jS^) and (^6|) as follows, 

S""^^^ = 2h" + 4^h' ^lAh + 2a^Spr (43) 



Here is given by Eq. (H) such that = - (a'/a) The integration function S^'^^ = ^g^^, (x^) can 
be further determined by substituting S"'^ in Eq. (^. Upon inserting these expressions in Eq. ( p^ ) we 
get the following inhomogeneous equation for S'^: 

5- " + 2^^^"/ - AS% = ^o"^^^, + ^0,^"'^ + h-$ + \6-,Ah - 3 ^ ( h-\p + V " ) 

- / [A ( /^"O,^ + /i^o " ) + f aV. ( + h^C )] 

(h" + 2^h' + a''6pr) . (45) 
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5. k=0 universes 

The unperturbed {hab = 0) metric satisfies (Misner 1973) 



3I /day_p^_p^^Q (46) 



such that the density pi equals pmo or Pmr some prescribed conformal time rj — r]Q. Using the Hubble 
constant 

(which is realistic in a matter-dominated ambience at r/ = 779), this equation can be written 



daV 8 / 2 



where 



2 „ Prpao jj 

p - -7. — 



is a constant. 

Hence, fixing the origin of ry, the solution has the form 

a=j,{^'- ^?) . (47) 

Normalizing the time coordinate by 

P 

we have 

« = ie - 1) (48) 

so that the Big Bang occurs at ^ = 1. The range of the time coordinate is determined by the ratio of the 
radiation density to matter density]^, thus in our Universe, ^ > 1. 

The evolution of the trace h, Eq. (^), is driven by the velocity potential [Eq. (|26|)]: 

h" + -h' + ^ {2E{x^) -h)^ ~^apr^'\' (49) 
a Ha 

where Eq. (17) was used for the matter density perturbation. We may get rid of the inhomogeneous term 
E{x^) by introducing the function 

f ^h- 2E{x'^). (50) 

Gauge transformations alter / as follows, 

f = f + Q^b. (51) 
With the new time variable ^, Eq. (49) takes the form 

f + - = -^T^Kij (52) 



^Small values, |g| < 1, are compatible with a big crunch. 
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where an overdot means d/d^ and 
Equation (E3) becomes 



3/4 



We now introduce a new, gauge-invariant potential u by writing 

u = A't' + /. 

Equation ( p^ and the time derivative of Eq. (|5^) then are, respectively. 

2? i / 6 I 16 W _ 16 



e-i ' ie-i) 



f 



0. 



(53) 

(54) 

(55a) 
(55b) 



Eq. (55a). Thus we get the fourth-order equation for u 



From (55b), we can express A/ and substitute it in the equation obtained by acting with the Laplacian on 

;ion for i 

JL / ^ I 2; d _ _6 16 \ dfji — n (^FA 



\d^2 ^ 



Au - 



Given a solution of the coupled equations (55a) and (55b), the radiation density perturbation can be 
computed from Eq. (^6|) as follows: 

So - J-o^/'^v - J^o^/^^ 



(57) 



6. Particular solutions 

[i) First we seek solutions with a vanishing gauge-invariant potential, u = 0. By Eq. ( |55b| ), then / is a 



harmonic function. A/ = 0, and Eq. (55a) becomes 

/ 



2C 



-/ 



6 



16 



This is the generalized Legendre equation for i/ = 2 and n — 4. Particular solutions are 

ri 1 



(58) 



and the associated Legendre function of the second kind 



QtiO = (1 - e 



,2 d^Q2 



(59) 



Inserting here the Legendre functions P2 = \ (3^^ — l) and Q2 — \P2 In — we get the particular 
solutions in the form 



/i 

/2 



48P4 



1 - 5^^ + 15g^ + 5 
5 



{e - 1)^ 



(60) 



48^2 
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Thus the solution of Eq. (|5q) is 



(61) 



where the combination functions Ai and A2 are harmonic, AAi — AA2 = 0, to yield the property A/ = 
as required. We may then perform a gauge transformation ( |5T| ) with a harmonic b such that ^2 = is set 
in the solution (|61 



(m) Our second set of particular solutions arises from the assumption that Au = 0. In this case we see 
from Eq. ( ^6|) that it is it, rather than /, that satisfies the generalized Legendre equation (|58|). Hence 



u^Ai {x") ui + A2 (x") U2 + C ix") i + D (x") 
where the two independent solutions are 



(62) 



u, EE 48//P|d6d6 = 6TI^'-^f'-|ln(e'-l) 



U2 



It follows from Au — that each of the integration functions Ai (x") , A2 (a;") , C (a;") and D (a;") is 
harmonic. 



A particular solution of Eq. (55a) is given by (Bronstein 1974) 

^(C) 



m) = ley [/2(o/i(c) - /i(o/2(c)] dc 

Since u is a harmonic function, so is /q. Using (|60| ), we obtain 

/o(0 = i(^{[-^(^^~3)^^i + i(3e^-6e^-25)A2 + (e4-2e2^5)C] 

X "(^2^l)ln|^-2ein(e2-l)" 
+ [-^ - + 125) + i (3^2 - 9) ^^2 + {2e - 6) ^C] 

X "(e2 + l)ln(e2-l)-2einf^" 
- (3 - - + m^4) + i (7 - 14^2 + 3^^) 
+2(2^' + + 5)^C + 20 (1 + e') i?} . 



(63) 



The solution of the inhomogeneous Eq. (55a) has the form 

/ = /o + Fi(a;")/i+F2(x")/2. 



(64) 



Equation ( 55b ) has yet to be satisfied: 
Inserting here (|6^), we get 



ii+ii-A/^O. 



Ai (a;")/i+A2 (x")/2 + ^A/ = 0. 



Hence 



2 - ^AF2 



where Fi and F2 are biharmonic functions, i.e., AAFi = AAF2 = 0. 
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(in) We may continue the process of generating new solutions by replacing next the harmonic condition 
on u with a biharmonic condition. Thus the function A A/ satisfies the homogeneous equation (^8|). The 
solution of this will define Au by the equation 



Au = -^AA/. 



(65) 



We take the Laplacian of both sides of Eq. (55a) and solve for A/, given the source term Am. This in 
turn yields the potential u by using Eq. ( ^5b| ). A comparison of ( |6l] ) and (|6j) reveals how the solution 
generating procedure essentially proceeds: the next solution in the sequence is generated by relaxing the 
harmonic condition on the coefficients of the given solution /. 



The general solution of Eq. ( p9| ) is 

h = MO + MhiO + ^2/2(0 + m^^)- 



(66) 



where Ai = Ai{x^) and A2 = A2{x^) are integration functions. The gauge transformations (|^ with a 
nonvanishing c" parameter alter h as follows: h = h + 20°^^. Thus, by a suitable gauge transformation, we 
arrange that E — 0, and still we can perform transformations with c'^ = 0. 

Finally, we get the trace- free part S"^ of the metric perturbation from Eq. (^8|). The solution will have 
the form 



J fi — O f 



/3 



where 5' is a spheroidal wave function (Flammer 1957, stratton 1935, Fisher 1937), a solution of the 
homogeneous equation 



(67) 



and S* 'fl is a particular solution of (pq). 



Let us consider the solutions which are given by C°° functions. Following White 1973 , we may then 
represent the amplitude Ai in terms of a C°° function B as follows, 



Ai = AS. 



(68) 



The treatment of the curl terms containing is a fairly straightforward task. There remain to be found 
the pure density perturbations with /i"" =S =0. For these perturbations, Eq. ( p8| ) simplifies somewhat. 



'0 ,/3p 



•^0/3/., 



\5-pAh 



3 



For solution (i), we have h = A_B/i and Spr = whence 

+ 1^^^", - M^A5", = [s^^^^ + s,^^^^ + h; - f AB- / ^ / Ada 



(69) 



We seek a particular solution in the form 



5 7^ = AS 



iJ^AS X(0 + ^o"/3 



(70) 



(71) 
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where the requirement of compatibiHty with the divergence equation ( p4[ ) yields 



(72) 



15 



3r-i2? +1 



and F {^) is a function to be determined. By substituting (|7l|) into (|70|), we get the ordinary differential 
equation 



4|_^Q _ „8,,2 



/3 F^-' /3 



(73) 



This has the solution 

F = 



- lokf^ [-63^' - 1249^' + (45^^ + 544^^ - 709) ln{e - 1) 
hl20ein 1^ - 120(e' - 1) ln(e - 1) ln(^ + 1) 



(74) 



where ci and C2 are integration constants. The terms in 5 proportional to C2 have the harmonic 

amplitude AS'^. These terms solve the homogeneous equation. They are a special case of gravitational 
waves described by spheroidal wave functions. This reflects on the ambiguous nature of the decomposition 
of the perturbations into wave and non-wave parts. 



Collecting the results, the tensor perturbation of solution (i) has the form 
where S af3 and /i are given in Eqs. (71) and (|60|), respectively. 



(75) 



The Sachs- Wolfe effect 



The temperature variation 5T of the cosmic background radiation can be computed ( Sachs and Wolfe 



1967) as follows, 



where rjji and rjE denote the time of reception and emission, respectively, and w is the affine length along 
the null geodesic of propagation with tangent 

^ = (-l,e") (77) 
dw 

such that e"ett = —1. We consider the contribution of the relatively increasing mode. Then /iq^ = and 
the second term under the integral in Eq. (|7^) vanishes. The term 5* ap has the amplitude B^ap- By using 
the relation 

y^a-j—dw = y aC^dw - y'd-q, (78) 
dw 

we get dipole anisotropy contributions with respective amplitudes B^pe^ and AB^^e^ and gravitational 
redshift terms. The trace part of hafi, unlike that of the pure dust, is time-dependent, and thus the 
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cancellation of the integrated terms in the (Sachs and Wolfe 1967) result does not occur here. Taken 
together and excluding dipole anisotropics, the contributions to the temperature variation sum up to 



ST 
T 



= TR — TE 



where tr and te are the values of the function 



2 AB 
675 IC+TF^ 
2 AS /qf2 



'675 ii-l)-' 
1 B 

'15 (e-i)'-' ' 

-8 (5^1^)3 [C2£. 



(3f 



9C^ 



13^+15) ln(e- 1) 
13^-15) In (e + 1) 
15C^ - 15) 
^(18^6 _ 220e^ - 885^' - 225)] 



(79) 



(80) 



at the respective events R of reception and E of emission. In addition to these contributions representing 
the original Sachs- Wolfe effect, we have the integrated term 



At = 



AB 



3 (C2_i). 



dw. 



(81) 



At large values of = /x^, the function r increases logarithmically. The integrated term At gives a negligible 
contribution at late times, but it must be taken into account in the vicinity of f = 1. This indicates that 
the temperature fluctuations in solution (?) can be larger than for a pure dust cosmology. 



8. Discussion of the results 

The evolution of density contrast of the incoherent matter, 5pml Pm mirrors the evolution of the trace 
h, as can be seen from Eq. (p^. For late times, that is in a matter-dominated era, one should reasonably 
expect this density contrast to be well-approximated by the Sachs- Wolfe scenario. However, for solution 
(z), the trace perturbation with the coefficient Ai behaves asymptotically as /i oc unlike the relatively 



growing mode of Sachs and Wolfe 1967. The perturbation with the coefficient A2 is /12 c>c 1/^ , precisely as 



the relatively decreasing mode of Sachs and Wolfe 1967 . The trace- free part of solution (i) at late times 
diverges also faster than in the absence of radiation: it goes like hi cx ^^In^. The asymptotic behavior of 
solution {ii) is similar: the coefficients of A2 and C tend to stationary values, and the coefRcient of D is 
asymptotically oc 1/^^. 

One may ask whether or not the perturbative cosmologies obtained here are favored by the exact, 
nonlinear evolutionary processes. This issue may be addressed by numerical methods, to be discussed 
elsewhere ( Czinner 2001] ). 
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